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INTRODUCTION 


Particles  bombarding  a  single  crystal  have  been  found  to 
exhibit  anomalous  behavior  when  incident  along  low  index  crystal 
directions.     Energy  loss  has  been  shown  to  be  unusually  small  in 
these  directions  and  some  of  the  projectiles  penetrate  to  very 
great  depths  in  the  crystal.     For  many  processes,  such  as  sputter- 
ing and  backscattering,  which  require  a  close  encounter  between 
projectiles  and  the  atoms  of  the  target  crystal,  minima  occur  when 
the  projectile  beam  is  aligned  close  to  a  low  index  direction. 
These  phenomena  Indicate  that  part  of  the  incident  beam  of  parti- 
cles have  a  lower  limit  imposed  on  the  distance  between  themselves 
and  the  target  atoms.     This  can  be  explained  if  the  atoms  of  the 
rows  in  the  low  index  directions  provide  a  set  of  correlated  colli- 
sions with  the  projectiles  which  steer  them  into  the  large  open 
spaces  or  channels  between  the  rows.     In  this  way  a  portion  of  a 
beam  of  particles  suffers  only  small  changes  of  direction  and  tra- 
vels directly  into  the  crystal  with  a  relatively  small  loss  of 
energy. 

A  theory  of  channeling  has  been  developed  by  Lindhard  (1), 
which  predicts  the  maximum  angle  a  channeled  particle  can  make 
with  a  low  index  crystal  direction.     The  theory  treats  the  problem 
in  two  approximations  relevant  to  high  and  low  particle  energies. 

The  high  energy  expression  for  critical  angle  has  been  shown 
to  agree  with  experiment  by  Davies  (2) .     The  low  energy  critical 
angle  crudely  agrees  with  the  experimental  results  of  Andreen  and 
Hines  (3).    However,  a  close  examination  of  this  data  showed  that 
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the  dependence  on  energy  may  be  better  explained  by  an  alternative 
expression  which  is  examined  in  this  thesis. 

Onderdelinden  (^)  has  been  able  to  predict  the  sputtering 
yield  of  single  crystals  with  projectiles  incident  near  the  low 
index  crystal  directions.     Sputtering  is  assiamed  to  be  caused  by 
the  energy  loss  of  the  unchanneled  fraction  of  the  beam  in  some 
surface  layer  whose  thickness  is  an  adjustable  parameter.     On  the 
basis  of  the  success  of  Onderdelinden ' s  treatment,  minima  in  sput- 
tering yields  are  direct  evidence  of  channeling. 

In  the  present  work,  single  crystals  of  silver,  aluminiim,  and 
gold  were  bombarded  with  Ions  ranging  in  atomic  number  from  1  to  18 
and  energy  from  10  to  100  kev.     The  sputtering  yield  was  measured 
as  a  function  of  the  angle  between  the  ion  beam  and  the   <110>  and 
<100)  crystal  directions.     The  angular  half  widths  at  half  minimum 
of  the  sputtering  yield  were  measured.     The  dependence  of  these 
half  widths  upon  energy,  and  atomic  numbers,  of  the  incident  ions 
and  crystal  atoms  was  compared  to  the  values  predicted  by  the  low 
energy  theory  of  channeling. 


THEORY 


Some  fraction  of  a  perfectly  colllraated  beam  of  energetic 
ions  incident  upon  a  perfect  single  crystal  along  a  low  index  di- 
rection has  small  impact  parameters  with  the  crystal  atoms  at  the 
surface.     This  fraction  of  the  beam  suffers  large  angle  scattering 
and  interacts  with  the  crystal  in  a  random  way.     The  rest  of  the 
beam  enters  the  relatively  open  spaces  or  channels  bounded  by 
atomic  rows.     The  ions  in  this  fraction  of  the  beam  will  oscillate 
in  the  channels  as  they  travel  into  the  crystal.     It  is  convenient 
to  think  of  these  as  the  aligned  and  random  portions  of  the  beam. 
If,  as  Ions  enter  the  crystal  at  some  small  angle  to  a  channel  di- 
rection, the  energy  in  the  motion  transverse  to  the  channel  direc- 
tion is  less  than  the  potential  barriers  of  the  channel  boundaries, 
these  ions  will  be  trapped  in  the  channel  as  they  travel  into  the 
crystal. 

For  purposes  of  studying  channeling,  the  atomic  rows  may  be 
approximated  by  continious  strings.     The  potential,  U(r) ,  of  such 
a  string  is  the  sum  of  the  potentials  of  all  the  atoms  in  the  row 
averaged  over  the  atomic  spacing  in  the  row. 


r  is  the  perpendicular  distance  from  a  field  point  to  the  row,  d 
is  the  distance  between  atoms  in  a  row,  z  is  the  distance  along 


ion  and  an  isolated  atom  of  the  crystal.     The  collision  between  an 
ion  and  a  string  is  assumed  to  be  Independent  of  the  neighboring 
strings. 
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the  row,  and  V  (r^  +  z^/2is  the  interaction  potential  between  the 


This  string  approximation  requires  that  the  ion  undergo  cor- 
related collisions  with  the  atoms  in  the  row  in  the  process  of 
being  steered  back  into  the  channel.     If  an  ion  did  not  follow  the 
trajectories  of  the  channeled  ions,  it  would  interact  with  an  in- 
dividual atom  in  the  string  for  a  time 

r. 


At  = 


min 


V  sin^ 

••^min  distance  of  closest  approach  of  channeled  ions  to  the 

string,  and  ^ is  the  angle  between  the  velocity  of  the  ion  and  the 
channel  direction  at  the  center  of  the  channel.    The  validity  of 
the  string  approximation  requires  that  the  distance  traveled  by  an 
ion  parallel  to  the  channel  in  this  time  be  greater  than  the  spacing 
of  the  atoms  in  the  row.    Since  the  velocity  parallel  to  the  channel 
direction  is  v  cos  ^  ,  this  condition  gives  , 


For  small  ^  this  becomes  V  sin^ 

'Tnin     ^  ^2  J 

The  motion  of  a  channeled  ion  can  be  projected  onto  a  plane 

normal  to  the  string  direction.     The  trajectory  of  the  ion  in  this 

plane  will  have  an  impact  parameter  I  with  the  string.     For  any 

particular  values  of      and  ion  energy  E,  the  value  of  rj^^^^  will  be 

the  smallest  when    I  =0.    This  situation  is  the  most  restrictive 

case  of  (2)  and  requires  that  the  kinetic  energy  in  the  normal 

plane,  E^,  becomes  zero  at  r^^^.     The  expression  for  the  critical 

angle  is  developed  for  the    I  =0  case. 

The  energy  balance  for  transverse  motion  Is  given  by 
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♦  U(r)  : 
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At  the  minimum  distance  of  approach  and  for  small  angles  this 
becomes  ,  2 


(3) 


Using  some  appropriate  form  of  string  potential  and  the  inequality 
(2)  to  eliminate  t^^^,  it  is  possible  to  obtain  an  estimate  for  a 
maximum  angle  at  which  channeling  can  occur.     Lindhard  (5)  has 
used  this  method  to  develop  expressions  for  the  critical  angle  in 
the  cases  of  high  and  low  energy  channeling. 

Lindhard 's  treatment  of  high  energy  channeling,  based  on  a 
Thomas-Fermi  atom-atom  potential,  uses 

2  In  Ca 


7  7 

U(r)  =       ^2  ^ 


where  the  screening  length  a  is    (  Q-  6853) 


and  a^  is  the 


Bohr  radius.     This  is  a  reasonable  approximation  to  the  string 
potential  for  r  <a.     C  is  a  constant  of  integration  which  appears 
in  a  general  expression  for  the  string  potential  given  by  Lindhard 
and  is  assigned  the  value  •>yj7 

Using  the  relations  (2)  and  (3),  r^^^^  is  eliminated  with  the 
result 
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'    ^     2  In  Ca 


>  sr. 


This  can  be  rewritten  as 

Ca 


dxj; 


exp 


2b 


>  1 


772 
where    b  =  ^1^2" 

E  , 
it  will  be  violated  first  by  the  exponential  as      is  Increased. 


If  the  inequality  holds  for  very  small 


Extremely  small  values  of  the  exponential  can  be  eliminated  by 


requiring 


The  inequality  will  hold  then  if 


ca 

The  first  of  these  gives  an  estimate  of  the  critical  angle. 


>  e. 


2b 


2  ZiZoC 


Ed 


(^) 


The  second  provides  limitations  on  the  energy  region  in  which  the 


treatment  is  valid. 


e. 
C 


> 


2  Zi Zge^ 
Ed 


Since  ^y=^  >  1 

or     E  >  e'=  2Z,Z,e^d 

a2 


this  becomes 
(5) 


Smaller  energies  result  in  larger  values  of  t^^^.     For  this 
reason  Lindhard  uses  an  approximation  to  the  string  potential  that 
is  valid  in  a  region  r  >  Ca, 


U(r)  =  l^ZzE^ 


Ca 


Using  (2)  and  (3)  to  eliminate  r 


min' 


the  relation 


(Car 


b  >  xj;' 


is  obtained.     This  can  be  written  as 


Ca  d'2 

"2 

C  a  ^ 

_  dd^'2 

\f2  d  \ 

[Car  z^z^ 


(6) 


Since  the  approximation  to  the  string  potential  was  based  on  r >  Ca, 
it  follows  that  rmin>Ga.     In  so  far  as  (2)  becomes  t^^^  =4^2'^' 

relation  \//  >         gives  a  limitation  on  the  energy  region  to  which 

^  C 
this  treatment  should  apply.     Since  r  1  ^     the  last  inequality 

a  /  1  2  V 

can  be  written  as  ^U/  > 

d  M 

becomes 


2Z.,Zge' 
Ed 


^  2 


.  With  the  expression  for  \b  this 
2ZiZ2e2d 


>  E 


E  <  E  . 


t-r 
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Llndhard  has  developed  these  two  cases  to  give  estimates  of 
the  behavior  of  the  critical  angle  for  channeling.     The  expressions 
for  ijj  ,  and  xjj   are  critical  angles  in  the  sense  that  particles  in 
the  aligned  portion  of  a  beam  make  angles  ^  with  the  channel  axis 
such  that  ^  <  c'V(    or  ^  <  c'^  .    The  constant  C,  which  Llndhard 
uses,  is  between  1  and  2. 

Both  expressions  for  the  critical  angle  predict  dependence  on 
the  energy,  atomic  numbers  of  the  incident  particles  and  target 
atoms,  and  lattice  spacing  of  the  crystal.    These  treatments  deal 
with  interactions  occuring  at  ranges  such  that  r  < a  in  the  high 
energy  case  and  r > Ca  in  the  low  energy  case.    It  Is  interesting 
to  consider  the  effects  of  interactions  occuring  in  the  neighbor- 
hood of  r=a. 

Moliere  (6)  has  given  an  approximation  to  the  Thomas-Fermi 
potential  with  the  form 


V(r)  =  ^1  ^2^' 
R 


0.1  exp 

.6R 
a 

+  0-55  exp 


.  1.2  R 
a 


♦  0.35  exp 


-0.3R 


This  approximation  has  been  used  by  Erginsoy  (7)  to  develop  an  ex- 
presslon  for  the  string  potential. 


U(r)  =  ^^1^2^^ 


0.1  Ko 


6r 
a 


*  0.55  K. 


1.2  r 


*  0.35K. 


0.3  r 


(7) 


The  functions  Kq  which  appear  in  this  expression  are  modified 
Bessel  functions  of  the  second  kind.     In  the  neighborhood  of  r=a, 
the  string  potential  can  be  reduced  to 


U(r)  =  BTTaZ^Zge- 


rd 


(8) 
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with  the  use  of  polynomial  approximations  to  Kq,  Abramowltz  and 
Stegun  (8).     B  is  a  constant  on  the  order  of  0.5.     Using  this  ex- 
pression for  the  string  potential  and  the  requirements  for  channel- 
ing (2)  and  (3),  an  expression  for  the  critical  angle  was  developed. 


BTTaZiZge^ 
This  can  be  rewritten  as 


> 


Ed2 


(9) 


For  ease  of  comparison  the  expressions  for  "y    and  XU  are 

1  '2 

listed  again. 

"21^12^2 


E  d 


Sa^Z^Zge^' 


Ed- 


(6) 


to 


Whenever  Z-j^  is  much  less  than  Z2,  the  screening  length  reduces 

a  -  0.6853 


With  this  approximation  the  dependence  of  critical  angle  on  Zi  or 
Z2  has  a  much  more  simple  form. 

The  predicted  dependences  on  E,  Z^,  and  Z2  are  summarized  in 
Table  I  for  l//  ,  ^2  ,  and  t//^  . 


9 


EXPLANATION  OF  PLATE  I 


Table  I.    A  svunmary  of  the  behavior  of  the 
critical  angle  predicted  by  ^i, 
xjjZf  and 
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PLATE  I 


TABLE  I 

Dependence  of  ^^'^  V^c'  °^  ^»  ^: 


Critical  Angle 

Zl 

Z2 

z 

1/2 

^2 

g.iA 

Zl 

1/12 

Zl 

2/9 

Z2 
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SPUTTERING 

The  developments  of  the  critical  angles  for  channeling  have 
all  used  the  basic  assumption  than  the  channeled  particles  ap- 
proach no  closer  than  some  minimum  distance  to  the  atoms  of  the 
target.     This  assumption  is  well  Justified  by  experiments  which 
measure  secondary  processes  which  depend  on  the  close  approach  of 
a  projectile.    Such  experiments  have  shown  minima  in  backscatter 
yield  and  yields  for  various  nuclear  processes  when  particles  bom- 
barding a  single  crystal  were  incident  in  a  channeling  direction. 
Similar  minima  have  been  observed  in  sputtering  yields. 

Onderdelinden  (9)  has  been  able  to  calculate  sputtering  yield 
as  a  function  of  l//  by  assuming  that  it  is  proportional  to  the 
energy  lost  by  the  nonaligned  portion  of  a  beam  of  particles  with- 
in a  certain  depth,  Xq,  of  the  target. 

These  calculations  were  first  performed  for  a  beam  of  parti- 
cles parallel  to  a  channeling  direction,  with  the  proportionality 
constant  and  x^  as  adjustable  parameters.     The  results  were  com- 
pared with  measured  values  of  the  sputtering  yield  at  various  ener- 
gies in  order  to  fix  the  parameters.     Using  these  parameters  and 
estimates  of  the  nonaligned  portion  of  the  beam  based  on  Lindhard's 
low  energy  critical  angle,  Onderdelinden  was  able  to  calculate  the 
dependence  of  the  sputtering  yield  on  the  angle  between  the  channel 
and  the  beam.     These  calculations  were  found  to  be  in  good  agree-  . 
ment  with  experimental  data. 

The  success  of  Onderdelinden' s  method  of  calculating  sputter- 
ing yield  indicates  that  sputtering  yield  is  a  good  measure  of 
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channeling.     In  the  course  of  his  development,  Onderlinden  found 

o 

that  Xq  was  on  the  order  of  100  A,  so  that  sputtering  measures 
channeling  at  the  surface  of  the  crystal.     Backscatter  studies  by 
Da vies  (10)  indicate  that  the  critical  angle  is  a  function  of  depth 
and  temperature  In  the  target.    However,  the  critical  angle  was 
shown  to  be  insensitive  to  temperature  at  the  surface.  Therefore, 
sputtering  measurements  provide  a  method  to  study  channeling  at  the 
surface  of  a  single  crystal  and  eliminate  problems  associated  with 
temperature. 

Plate  II  illustrates  the  energy  and  temperature  dependence  of 
sputtering  yield,  Almen  and  Bruce  (11)  found  that  over  the  energy 
range  studied,  the  sputtering  yield  was  fairly  large.  This  allows 
measurements  to  be  made  with  low  beam  currents  to  minimize  crystal 
damage.  The  stability  of  the  sputtering  yield  observed  by  Nelson 
(12)  in  the  neighborhood  of  ambient  temperatures  is  another  factor 
in  the  desirability  of  sputtering  measurements  from  the  experimen- 
tal point  of  view. 
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EXPLANATION  OF  PLATE  II 


Fig.  1,     The  sputtering  yield  of  Ag  bombarded 
by  Ar,  Ne,  and  N  at  various  energies. 

Fig.  2.     The  behavior  of  the  sputtering  yield 
of  Ag  with  temperature. 


PLATE  II 
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EXPERIMENTAL  METHOD 

Sputtering  measurements  were  made  with  Au,  Al  and  Ag  single 
crystals.     These  crystals  were  produced  by  Semi-elements,  Inc., 
and  were  grown  to  a  size  of  about  0.5"  x  0.5"  x  0.25"  with  the 
<111>  direction  normal  to  the  surface.    All  of  these  crystals  have 
fee  structure  with  lattice  spacings  of  ^.0^^,  ^.07^,  and  ^.08^  for 
Al,  Au,  and  Ag  respectively. 

The  crystals  were  rotated  about  an  axis  parallel  to  the  <111> 
direction  as  shown  in  Plate  III.     This  rotation  was  monitored  and 
recorded  as  the  angle  cj)  .     The        axis  was  set  at  an  angle,  ©  ,  to 
the  ion  beam  equal  to  the  angle  between  the  <111>  direction  and  the 
channel  direction  to  be  studied.     The  cj)  rotation  was  continuously 
driven  by  a  motor  so  that  the  desired  channel  direction  was  ro- 
tated through  the  ion  beam. 

The  charge  given  off  by  the  crystal  was  collected  with  a 
large  plate  subtending  a  solid  angle  of  about  27r  steradians  at 
the  crystal  surface.     The  current  between  the  collector  and  ground 
was  measured  by  a  Keithley  Picoammeter.     The  collector  plate  could 
be  biased  with  respect  to  the  crystal.    With  the  ion  beam  striking 
the  crystal  in  a  random  direction,  the  measured  current  was  record- 
ed for  various  biasing  voltages.     This  data  was  plotted  as  current 

yield,  measured  current  divided  by  incident  ion  current,  versus 

4  + 

biasing  voltage,     Plate  IV  shows  the  results  for  a  beam  of  He 
ions  bombarding  a  gold  crystal.     The  positive  yield  quickly  stabi- 
lized to  a  value  of  R"*"  on  the  order  of  1,  near  a  biasing  voltage 
of  -20  volts,    r"*"  could  be  composed  of  backscattered  positive  ions, 


EXPLANATION  OP  PLATE  III 


A  schematic  representation  of  the  orientation 
of  the  crystal  in  the  target  chamber  with  re- 
spect to  the  ion  beam  and  the  collector. 
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PLATE  III 
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secondary  electrons  ejected  from  the  collector  plate  by  the  ener- 
getic backscattered  particles,  or  positively  charged  sputtered 
target  atoms.     However,  estimates  of  the  backscatter  yield  in  this 
energy  region  are  on  the  order  of  10"^.     Therefore,  backscattered 
Ions  and  their  secondary  electrons  can  be  Ignored  in  R^.    All  sput- 
tering measurements  were  taken  at  a  collector  bias  of  about  -70 
volts. 

Sputtering  yield,  R"*",  was  recorded  as  a  function  ofcj),  on  a 
Brush  strip  chart  recorder,  and  minima  were  observed  to  be  centered 
about  cj>  values  that  correspond  to  low  index  crystal  directions. 
This  data  was  taken  for  the   <L00>  and  <110>  directions  in  all  three 
crystals.    Beam  currents  and  sputtering  currents  were  on  the  order 
of  10"^°  amperes.     The  crystals  were  bombarded  with  ions  from  H"*" 
to  Ar"*"  with  energies  from  10  to  100  kev.     The  energy  of  the  ion 
beam  was  uncertain  by  2^. 

Examples  of  the  sputtering  yield  minima  are  shown  in  Plate  V. 
Three  distinct  minima  were  observed  for  each  revolution  of  the 
crystal.     This  is  caused  by  threefold  symmetry  about  the  <111> 
direction.     Ninety  degree  increments  of      were  marked  on  the  strip 
chart  recorder.     The  linear  half  width  at  half  minimum  of  the 
sputtering  minima  were  measured  and  converted  to  an  angular  measure 
of  ^  . 

According  to  Onderdelinden' s  treatment  of  sputtering,  the 
half  minimum  of  the  sputtering  yield  corresponds  to  a  situation  in 
which  the  aligned  fraction  of  the  beam  is  one  half  of  that  value 
it  would  have  for  ^=0.     The  sputtering  yield  apptoaches  the  base 
line  for  random  direction  sputtering  as  X^f  approaches  the  critical 


EXPLANATION  OF  PLATE  V 


Three  examples  of  the  recorder  output  with  90° 
incriments  of  (j)  indicated.     The  large  minima 
represent  channeling  in  the  direction  consid- 
ered in  each  case.     The  smaller  minima  corre- 
spond to  other  channeling  directions, 

Pig.  3.     Collected  current  for  45  kev  H"*"  ions 

incident  upon  the  <100>  direction  in  Ag. 

Fig,  4,    Collected  current  for  25  kev  Ne"^  ions 
incident  upon  the  <100>  direction  of  Ag, 

Fig,  5.     Collected  current  for  100  kev  N"*"  ions 
Incident  upon  the  <110> direction  of  Ag. 


22 


23 


angle.     In  addition  the  angular  divergence  and  small  energy  spread 
of  the  ion  beam  make  the  intersection  of  the  minimum  with  the  base 
line  an  unreliable  measure  of  the  critical  angle.     The  width  of 
the  sputtering  yield  minimum  at  half  minimum  is  distinct  and  rel- 
atively insensitive  to  the  angular  divergence  and  energy  spread  of 
the  beam.     For  these  reasons  the  critical  angle  was  not  measured 
directly,  but  the  half  width  at  half  minimum  of  the  sputter- 

ing yield,  was  used  to  determine  its  dependence  on  E,       ,  and  Z2 , 

It  was  observed  that  the  three  minima  in  a  single  rotation 
did  not  have  the  same  shape.     This  indicated  that  the      axis  was 
not  truly  parallel  to  the   <111>  direction  In  the  crystal.    As  the 
crystal  was  rotated,  each  of  the  three  equivalent  channeling  di- 
rections made  slightly  different  angles  with  the  ion  beam.  This 
problem  was  attributed  to  errors  in  mounting  the  crystals  and  in- 
accuracies In  the  crystal  surface.     It  was  necessary  to  observe 
the  behavior  of  one  of  the  three  minima  for  various  ©settings, 
A  value  of  ©  was  selected  for  this  minimum  between  extremes  of 
misalignment  with  consideration  for  depth,  width,  and  symmetry  of 
the  minimum. 

Several  measurements  of  the  sputtering  yield  were  taken  for 
each  selected  channel  at  a  particular  energy.     The  measured  half 
widths  were  averaged  to  obtain  mean  values  which  were  foiind  to  be 
reproducible  to  within  3^.     The  angular  widths  of  the  sputtering 
minima  in  terms  of  the'  angle  ^  between  the  ion  beam  and  the  chan- 
nel direction  were  calculated  from  the  relation  ^1  sin  ®  , 
This  relation  Is  easily  derived  from  the  geometry  of  the  crystal 
mounting,                                           ■  -     ■       -    ? -  ■ 
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APPARATUS 

The  main  elements  of  the  apparatus  are  shown  in  Plate  VI. 
The  desired  ions  were  produced  from  their  elemental  gaseous  state 
in  the  ion  source.     These  ions  were  drawn  into  the  accelerating 
tube  by  an  extraction  electrode  and  accelerated  to  a  predetermined 
energy.     The  ion  beam  was  defined  by  a  1  mm  slit  S^.    Since  the 
energy  of  the  ions  was  known,  the  bending  magnet  was  used  to  select 
the  desired  mass  by  assuming  charge  states  of  +1. 

The  ion  beam  was  collimated  by  two  1.5  mm  slits,  S2  and  S-^, 
separated  by  a  distance  of  one  meter.     The  angular  divergence  of 
the  beam  was  less  than  ±0.09°.     The  pressure  above       was  about 
10"^  torr,  and  between  S2  and  S-^  it  was  about  10~^  torr.     The  vac- 

o 

uum  in  the  target  chamber  was  maintained  at  about  10"    torr  by  the 
two  Ultek  Differential  Sputtering  ion  pvunps.    A  Faraday  cup  behind 
S-^  could  be  moved  into  the  beam  to  measure  the  incident  ion  current. 
The  cj)  rotation  was  monitored  by  the  use  of  four  electrical  contacts 
which  sent  a  voltage  pulse  to  the  recorder  at  each  90°  increment. 
An  adjustable  pot  well  insulated  from  ground  was  used  to  bias  the 
collector.     The  Keithley  Model  4l6  Picoammeter  has  a  rise  time  of 
less  than  0.04  sec.  in  the  current  range  used,  and  is  equipped 
with  a  3  volt  recorder  output.     This  was  used  to  drive  the  Brush 
strip  chart  recorder,  so  that  a  continuous  record  of  the  sputter- 
ing current  with  each  90°  increment  of  d)  denoted  was  obtained. 


EXPLANATION  OF  PLATE  VI 


The  main  features  of  the  apparatus  used  to 
make  the  ion  beam  and  measure  the  sputtering 
yield. 
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PLATE  VI 


Recorder 


Bending 
Magnet 


Target  Chamber 
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RESULTS 

The  observed  dependence  of         on  energy  for  various  pro- 
's' 

jectiles  in  Ag  is  displayed  in  Plates  VII  and  VIII.     If  the  de- 
pendence is  of  the  form         ocE"^,  with  p  a  constant,  a  log-log 
plot  of  ^1  versus  E  will  be  a  straight  line  of  slope  -p.     It  was 
found  that  straight  lines  were  a  good  fit  to  this  data  in  all 
cases.     The  slopes  are  listed  in  Table  II.     The  energy  E'  which 
divides  Lindhard's  development  into  high  and  low  energy  regions  is 
also  listed  in  Table  II.    A  3%  uncertainty  in  Xpx  makes  the  slopes 
uncertain  by  *  0.03. 

TABLE  II 

Measured  slopes  of  log-log  plots  of 
versus  E  and  calculated  values  of  E' . 

Projectile  Slope  E' (kev) 
100  direction 

h"*"  -0.38  3^6 

He"*"                               .        -0.38  885 

N*  -0.29  2890 

Ne"^  -0.30  •  ^360 
110  direction 

He"^              '        \  -0>2  r  -      ^    .  626 

N*  -0.33  20^3 

Ne"^  -0.32  ;v  3081 

Ar  -0.2^  4360 
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EXPLANATION  OF  PLATE  IX 

The  dependence  of  on  In  Ag,  The  data 
is  displayed  on  log-log  plots, 

O     Denotes  25  kev 

A     Denotes  50  kev 

□    Denotes  75  kev 
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PLATE  IX 
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EXPLANATION  OF  PLATE  X 

The  dependence  of         on  Z2  in  the  <110>  direc- 
tion of  Ag.     The  data  is  displayed  on  log-log 
plots. 

Fig.  8.  He"*"  incident  ions. 
Fig.  9.     H"*"    incident  ions. 

O     Denotes  25  kev 

A    Denotes  50  kev 

□     Denotes  75  kev 
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PLATE  X 
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EXPLANATION  OP  PLATE  XI 

The  dependence  of         on       in  the   <100>  dlrec 
tion  of  Ag.    The  data  is  displayed  on  log-log 
plots. 

Pig.  10.  He"^  incident  ions. 
Pig.  11.    H"*"    incident  ions. 

Denotes  25  kev 

Denotes  50  kev 

Denotes  75  kev 
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The  energy  region  studied  was  in  all  cases  below  E* ,  and  the 
observed  slopes  fall  in  a  range  between  -0.24  and  -0.44.  The 
steepest  slopes  were  observed  for  H"*"  and  He"*^  which  also  have  the 
lowest  values  of  E' , 

The  assumption  that  Z-^  is  much  less  than  Z2  made  it  possible 
to  describe  the  dependence  of  ^c  ^P°^  '^2 

a  simple  power  law.     This  predicted  behavior  was  summarized  in 
Table  I.     The  slopes  of  straight  lines  fitted  to  log-log  plots  of 

vs  Z-^  and  ^1  vs  Z2  were  measured  to  study  the  dependence  of 
the  critical  angle  on  and 

The  measured  values  of         are  plotted  against  Z^  in  Plate  IX 
for  various  ion  energies  in  Ag,    Straight  lines  provided  a  good 
fit  to  this  data,  and  the  slopes  are  listed  in  Table  III.  These 
slopes  are  uncertain  by  ±0.02  for  the   <110>  data  and  ±0.03  i'or 
the  <100>  data. 

TABLE  III 
Measured  slopes  of  ^1  versus  Z-t 


•  2  ^ 

Energy  (kev)  Slope 

<100)  direction 

25  0.39 

50  '  0.39 

75  0.41 

ai0>  direction 

25  0.37 

50  "  0.42 

75  0.43 
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The  Z2  dependence  of  the  critical  angle  was  studied  with  H**" 
and  He"*"  Ions  bombarding  Au,  Ag,  and  Al.     Log-log  plots  of  l//^^ 
versus       appear  in  Plate  X  and  Plate  XI.     The  slopes  of  straight 
lines  fitted  to  the  three  points  obtained  for  each  energy  are 
listed  in  Plate  XII,     The  slopes  of  these  lines  were  uncertain  by 
=t0.03. 

The  spacing  of  atoms  in  the   <110>  row,  d^^Q,  is  equal  to  (2)^ 


o 


f  the  spacing  in  the  <100>  row  d^QQ.     The  ratio  of  the  critical 
angles         /^<^qq>  "^^^  ^1'  ^2 .  and  E  fixed  is  (2)^/^  or  about 
1.26.     These  ratios  were  calculated  for  h"*",  He"*",  n"*",  and  Ne"*"  on  Ag 
for  energies  of  25,  50,  and  75  kev.     The  results  are  svunmarlzed  in 
Table  IV,  and  the  average  value  was  found  to  be  1.34. 

TABLE  IV  , 
Measured  values  of  the  ratio        /  t<100> 

Energy  (kev)  Projectilei      h"*"         He"*"  Ne"*" 

25  •  1.32      1.39      1.39  1>0 

50  1.22      1.3^      1.36  1.39 

75  1.23      1.31      1.3^+  1.38 

In  order  to  compare  the  three  expressions  for  the  critical 
angle  with  the  observed  behavior  of  the  predicted  dependences 

of   V'l.    %.  and         on  Zj^,  Z2  and  E  are  siunmarized  In  Table  VI 
along  with  the  range  of  observed  dependences. 


EXPLANATION  OP  PLATE  XII 


Table  V,    The  slopes  of  the         vs  Z2  plots. 

These  slopes  were  measured  from  the 
graphs  in  Plate  IX  and  Plate  X. 


PLATE  XII 


TABLE  V 

Measured  value  of  the  slope  versus 


Energy  (kev)  Slope 

0.10>  direction  -    "  ^ 

-      '  ' 

25  0.19 

50  0.16 

75                     >  0.15 

+ 

He 

25  0.15 

50  ^  ;              .  0.16 

75  0.17 

<100>  direction 

*25       .  0.07 

50  :    .'  0.09 

75  0.10 


25  0.08 

50  0.10 
75  0.11 
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CONCLUSIONS 


The  experimental  results  show  that  in  the  energy  range  stud- 
led,  the  critical  angle  for  channeling  can  be  described  by  a  re- 
lationship of  the  form   ^^^^o:  E"^  for  various  values  of  p.     It  was 
found  that  p  tended  to  decrease  as  E' ,  the  dividing  point  between 
Lindhard's  high  and  low  energy  regions,  increased. 

The  data  summarized  in  Table  VI  show  that  for  the  range  of 
ions  studied,   W    best  describes  the  energy  dependence  of  the  crit- 

leal  angle.     In  the  extreme  cases,  the  energy  dependence  was  found 

-0  2i^  -0.^^ 
to  be  E     '      and  E  .     These  values  bracket  the  behavior  pre- 

dicted by  \b  ,  but  are  outside  the  estimated  loncertainty  ±0.03. 
■  c 

These  results  indicate  that   l//    gives  a  reasonable  approximation 
to  the  energy  dependence  of  the  critical  angle  for  the  ions  stud- 
led  in  the  energy  range  10  to  100  kev. 

The        dependence  of  the  critical  angle  was  found  to  be  pro- 
portional to  ^i^^ .     The  observed  values  of  -p^  were  all  between 
0.37  and  0.^3.     This  behavior  is  about  midway  between  0.33  pre- 
dicted by    ^  and  0.50  predicted  by        .     However,  both  of  these 
predicted  values  are  outside  the  estimated  uncertainty  which  was 
±0.03. 

Only  crude  estimates  of  the  Z2  dependence  of  the  critical 
angle  were  possible.     Straight  lines  fitted  to  the  three  points 
plotted  for  each  case  showed  a  definite  difference  in  slope  for 
the   <100>  and   0.10)  directions.     Assuming  a  dependence  of  the  form 
Z2^^,  the  observed  values  of  pg  in  the   <3.00>  direction  were  be- 
tween 0.07  and  0.11.     These  agree  with  the  dependence  predicted 


by  0-08,  within  the  estimated  uncertainty  ±0.03.     "j^e  range 

of  P2  in  the  <110>  direction  was  0.15  to  0.19.  These  values  fall 
closer  to  the  prediction  of  which  is  0.22.  However,  only  the 
largest  is  within  the  estimated  uncertainty. 

It  was  found  that  none  of  the  expressions  for  the  critical 
angle  could  consistantly  predict  the  dependence  on  the  atomic 
numbers  of  the  incident  ions  and  crystal  atoms. 
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The  angular  widths  of  minima  in  sputtering  yield  were  used 
to  study  the  critical  angle  for  channeling  in  the  <100>  and  <110> 
directions  in  fee  crystals.     Silver,  aluminum,  and  gold  single 
crystals  were  bombarded  with  ions  ranging  from  H'*'  to  Ar"*^  at  ener- 
gies between  10  and  100  kev.     The  dependence  of  the  critical  angle 
on  incident  ion  energy  E  and  the  atomic  niimbers  of  the  incident 
ions  and  the  crystal  atoms  was  measured  and  compared  with  theoret- 
ical predictions. 

Three  expressions  for  the  critical  angle  were  given.     Two  of 
these,  and  were  developed  by  Lindhard  for  cases  of  high 

and  low  energy  channeling.  The  third  was  based  on  the  Moliere 
atom-atom  potential. 

The  crystals  were  rotated  about  the   <111>  direction  which  was 
normal  to  the  crystal  surface.     The  <111>  direction  was  set  at  an 
angle  to  the  ion  beam  equal  to  the  angle  between  the  <111>  direc- 
tion and  the  channeling  direction.     The  charge  liberated  from  the 
crystal  under  ion  bombardment  was  used  to  measure  the  sputtering 
yield  as  a  function  of  the  angle  between  the  ion  beam  and  a  chan- 
neling direction.     The  angular  half  widths  at  half  minimum  of  the 
minima  in  the  sputtering  yield  were  plotted  on  log-log  graphs  to 
measure  a  power  law  dependence  of  the  critical  angle  on  E,  Z^,  and 
Z2.     Straight  line  fits  to  these  plots  showed  that  a  power  law  was 
a  good  description  of  the  E  and       dependence.     The  three  points 
of  each  Z2  dependence  plot  were  in  reasonable  agreement  with  a 
straight  line. 

The  observed  powers  in  the  energy  dependence  fell  between 
-0.2^  ±  0.03  and  -OA'^  *  0,03.     These  were  best  predicted  by  l/z^. 
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with  a  dependence  of  E"^'^^.     The       dependence  in  the  <100> 
channel  direction  was  in  good  agreement  with  l/Zg.     In  the  <110> 
direction  the        dependence  was  found  to  be  midway  between  the 
behavior  predicted  by      ^  and  y^^.     Thus  the  dependence  of  criti- 
cal angle  on       and       cannot  be  consistently  described  by  any  of 
the  three  theoretical  expressions  considered. 


